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Abstract 

We discuss the new model expression Q;an(Q^) recently obtained for the QCD 
running coupling with a regular ghost-free behavior in the "low Q^" region. Being 
deduced from the standard "asymptotic-freedom" expression by imposing the Q^- 
analyticity - without any adjustable parameters - it obeys nice features: (i) The 
universal limiting value Q:an(0) = 47r//?o — 1.4 expressed only via group symmetry 
factors and independent of experimental estimates on the running coupling as{Q'^) 
(of QCD scale parameter A). This value turns out to be stable with respect to higher 
order corrections; (ii) Stability of IR behavior with respect to higher-loop effects; 
(iii) Coherence between the experimental aan(-^T) value and integral information 
on IR OLsiQ"^) behavior as extracted from jet physics data. 



The ghost-pole problem in the behavior of a running coupling, being an obvious prop- 
erty of the geometrical progression, spoils a physical discussion of the RG-summed per- 
turbative QCD results in the infrared (IR) region. To avoid it, one uses some artificial 
constructions like the "freezing of the coupling" hypothesis. 

Here, we are going to revive an old idea of combining the RG summation with an- 
alyticity in the variable. It was successfully used in the late 50 's for examining the 
QED ghost-pole issue [0, . Quite recently, it has been proposed for applying to the QCD 
case P]. 

The QED effective coupling q;(Q^) being proportional to the transverse dressed-photon 
propagator amplitude according to general principles of local QFT satisfies the Kallen- 
Lehmann spectral representation and, therefore, is an analytic function in the cut complex 

plane. 

The "analytization procedure" elaborated in papers |I|, |^ consists of three steps: 

(I) To find an explicit expression for aRciQ"^) in the Euclidean region > by 
standard RG improvement of a perturbative input. 

(II) To perform the straightforward analytical continuation of this expression into the 
Minkowskian region ReQ^ < 0, \m.Q^ = — e. To calculate its imaginary part and to define 
the spectral density by PRcicr, a) = Ima/jcl^cr — ie, a). 

(III) Using the spectral representation [see Eq. (Q) below] with p^c in the integrand 
to define an "analytically-improved" running coupling Q;an(<5^) in the Euclidean region. 

Being applied to a{Q^) in the one-loop ultra-violet (UV) QED case, this procedure 
produced 0] an explicit expression with the following properties: 

(a) it has no ghost pole, 

(b) as a function of a at the point a = it possesses an essential singularity 
~ exp(— Stt/q;), 

(c) in the vicinity of this singularity for real and positive a it admits a power expansion 
that coincides with the perturbation one (used as an input), 

(d) it has the finite UV limit a(oo, a) = Sir that does not depend on the experimental 
value a ~ 1/137. 

The same procedure being applied to the two-loop QED case yielded a more com- 
plicated expression with the same essential features. 

In the QCD case, to apply this technique to the strong running coupling, one has to 
make two reservations. 

First, as far as here asiQ"^) is defined via a product of propagators and a vertex 
function, there is a question about validity of the spectral representation. Happily, this 
point has been discussed in paper As a result, one can use the Kallen-Lehmann 
analyticity here, as well. 

Second, in QCD, the running of coupling is, generally, connected with the running of 
gauge. For simplicity, we assume that the MS scheme is used (or the MOM scheme in the 
transverse gauge) when ag is not infiuenced by the running of gauge. 

To construct an analytic effective coupling in the QCD case, we start with the leading- 
logs expression 

' 1 + aPo HQVf^') Po ln(gVA2) ^ ' 

with a{Q^) = ds{Q^)/4:n and Pq = 11 — (2/3)n/, the one-loop coefficient, and with the 
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spectral representation 

1 r°° , p(cr,a) 

aan(Q') = -/ da ■ ■ (2) 

vr Jo a + Q'^ — le 

According to step (II) of the outlined procedure, we define the spectral function in the 
one-loop approximation 

(1) / X a^PoTT a 

(1 + afjQL) + [ajjQ'KY 

Note that the RG invariance of d^n defined via Eq.(0) is provided by the scaling property 
of the spectral function 

p{a/n'^,a) = p(ln(T/A^), A^ = p^ exp(-l/a/?o) • (4) 

Substituting into Eq. (Q) we get [§] 



-am') = ^ 



1 A 

+ 



2 



lnQVA2 A2-g 



(5) 



where we used the QCD scale parameter defined as in Eq.(^. However, to identify a 
with a(p^), the running coupling value at = p^, we have to change this definition for 
A^ = p^ exp[— (/)(a(p^)/3o)], where the function (j){z) satisfies the equation z = l/4>{z) — 
l/[exp(0(^)) - 1]. 

It is clear that the "analytic" coupling constant, Eq. has no ghost pole at = A^, 
and its IR limiting value a^p{0) = 'i'n / [3q depends only on group factors. Numerically, for 
n/ = 3, we have ain^O) = 4vr/9 ~ 1.398. 

Usually, we are accustomed to the idea that theory supplies us with a set of possible 
curves for as{Q) and one has to fix the "physical one" by comparing with experiment. 
Here, Eq. describes a family of possible curves for aan(<5^) forming a bundle with the 
same common limit at = as it is shown in Fig. 1. 

Another important virtue of Eq. (|^) is that the analytic behavior in the IR domain is 
provided by a nonperturbative contribution ~ exp(— l/a/3o)- 

To analyze the two-loop case, let us start with aP'\Q'^) written down in the form 

where hi = /3i/(/3o)^ and Pi = 102 — 38/3r;,j is the two-loop coefficient. This expression 
corresponds to the result of exact integration of the two-loop differential RG equation 
resolved by an iteration and generates a popular two-loop formula with \nl/P term. 
For the spectral density, we get 



with 



/?opgk^) = ^,^^y^7^, ^ = ln^ (7) 



R{L) = L + bi\n 
I{L) = 71 + 6iarccos 



1 + ^ + ' 



\lV^ViJ ^Kbi 
bi + L 

(6i + L)' + 7r2 
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Figure 1: The bundle of analytic solutions. The curves (a) and (b) are the 1-loop an- 
alytic as for A =200 MeV and 4OO MeV, respectively; the curves (c) and (d) show the 
corresponding perturbative result. 



Now, to obtain a^nlQ^)' has to substitute Eq. (|^ into the r.h.s. of Eq. (H). 

However, the integral expression thus obtained is too comphcated for presenting in an 

—(2) 

exphcit form as the integration resuh differs from the used input aj^Q not only by the pole 
term "subtracting" the ghost pole (as in the one-loop case), but also by an integral along 
the unphysical cut < < A^exp(— 61) "born" by the log-of-log dependence. For a 
quantitative discussion we have to use numerical calculation. 

Nevertheless, for a particular value at = we can make two important statements. 
First, the IR limiting coupling value a;an(0) , generally, does not depend on the scale 
parameter A. This is a consequence of RG invariance (compare with Ref. 0) and in our 
case follows directly from Eq.(^). Second, the IR limiting coupling value is defined by 
the one-loop approximation, that is in the two-loop case coincides with the one-loop case: 

o^^M = «i!.^(o) =«,„(o). 

To obtain a simple proof for the two-loop case, it is convenient to express the difference 
•^in (0) ~ "^in (0) the imaginary part of the integral 

Jc z z + m(l + z) 

with the contour C defined by —00 < Kez < +00, Imz = —-k /hi. As far as the integrand 
is an analytic function in the half-plane below the contour C, we conclude that A = 0. 
Note also that the universality of aan(O) follows directly from the procedure of constructing 
the analytic coupling. Indeed, our two-loop input Eq. (||) has the ghost pole at = Is? 
and the unphysical cut mentioned above. The analytization procedure removes these 
parasitic singularities by two compensation terms. The term that removes the pole gives 
the contribution Aa(0)poie = l/(2/?o)- The second term that compensates the cut can be 
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Figure 2: "Higher loop stability" of the analytic solution. The normalization point is on 
the T lepton scale: as{M^) =0.34- Practically, the curve a^{Q'^) coincides with a'^^(Q'^) 
and we do not plot it in the figure. 



which equals l/(2/3o) and we again obtain the universal value of aan(O) = l/Po- Thus, in 
contrast to perturbation theory, where the many-loop corrections change the IR behavior 
of the running coupling significantly, our analytic coupling has a stable IR limit. The 
analytization procedure removing all unphysical singularities leaves us with the physical 
cut which is mainly described by the one-loop contribution. 

Note also that the universality of the value of aan(O) is not simply a matter of approx- 
imate resolution, Eq. (^, of the exact RG solution - it has a deeper ground. This fact 
can be established in a more general context, e.g., by considering the analytic properties, 
given by the analytization prescription, in the complex 1/a-plane. The details of this 
reasoning are rather lengthy and will be published elsewhere. 

Thus, the aan(O) value, due to the RG invariance, is independent of A and, due to 
the analytic properties, independent of higher corrections. This means that the causality 
(=analyticity) property brings the feature of the universality. 

Here, we mean also that the whole shape of the aa.n{Q^) evolution turns out to be 
reasonably stable with respect to higher corrections. The point is that the universality 
of Q;an(0) practically gives rise to stability of the a'^^iQ'^) behavior with respect to higher 
correction in the whole IR region. On the other hand, this stability in the UV domain is 
a reflection of the property of asymptotic freedom. As a result, our analytic model obeys 
approximate "higher loops stability" in the whole Euclidean region. Numerical calculation 
(performed in the MS scheme for one-, two-, and three-loop cases with = 3) reveals 
that a^^iQ^) differs from aW(<5^) within the 10% interval and a^XQ^) ^om a^JiQ^) 
within the 1% limit. This fact is demonstrated in Fig. 2 



expressed via the discontinuity of function (^) on this cut and presented as 




(10) 
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It seems that the IR stabihty is an intrinsic feature of a non-analytic (in a^) con- 
tribution. To illustrate this thesis, consider a recent IR modification, asiQ"^), for the 
QCD running coupling by Badalian and Simonov 0. These authors have studied a 
non-perturbative contribution to the QCD running coupling on the basis of the general 
background formalism using nonperturbative background correlators as a dynamic input. 
They came to the conclusion that these effects can practically be described by introduc- 
ing an effective gluonic mass Mb = V27Ta defined by string tension a into all "gluonic 
logarithms": InQ^ — >■ ln(Q^ + M^). This means that essentially slows down its 

evolution (i.e., freezes) around Q ^ Mb- Their numerical estimate gave Mb — 1.5 GeV. 
In practice, this yields the difference between one- and two-loop results in the IR region 
of 15% order. 

In our calculations, we used n/ = 3 as an average quark number. This seems to be 
reasonable in the low energy region Q < 2 GeV ~ rric. For a more realistic description of 
the as{Q) evolution in the whole Euclidean domain, one should take into account quark 
thresholds. To this end, one usually applies a matching procedure, changing abruptly the 
number of active quarks at an "effective" threshold Q = with some matching param- 
eter ^ — 1. Evidently, any procedure of that sort violates the analytic properties. On 
the other hand, these properties could be preserved by the "smooth matching" algorithm 
devised on the base of an explicitly mass-dependent RG formalism ascending to Bo- 
goliubov. This algorithm has been recently used for the precise analysis of the otsiQ) 
evolution in the 3 GeV < Q < 100 GeV interval. 

In the present context it could be applied also in the Q < A domain. However, this 
would change the a^niQ"^) behavior in the "very low Q" region only slightly as far as the 
limiting Q;an(0) value depends on Po{n*jr) with the effective quark number n*j- related to 
the ghost-pole position A. As it is generally accepted on the base of DIS data, in the MS 
scheme A„^=3(2 — loop) ^ 250 ^ 350 MeV which is quite above of the strange quark mass. 
This means that the use of value rij = 3 is justified. 

Analytic properties of the running coupling are important from the point of view of 
phynomenological applications, for example, for the description of the inclusive decay of 
the r lepton. To this end, one usually transforms the initial expression for the Rr ratio 
to the integral form in the complex g^-plane (see, e.g., [^]). This transformation based 
on analytic properties mentioned above which are violated in the standard perturbative 
consideration and maintained within our method. Note also that referring to "low-Q" 
data, like those of r-lepton decay, one should distinguish between QCD scale A'^^ in 
the usual RG solution taken in the MS scheme and A'^" corresponding to our analytic 
expression. For instance, in the one-loop case, to the as(M^) = 0.34 there correspond 
A|° = 280 MeV as compared with Af^ = 230 MeV. 

The idea that the QCD running coupling can be frozen or finite at small momenta 
has been considered in many papers (see, e.g., the discussion in \TU\). There seems to be 



experimental evidence in favor of this behavior of the QCD coupling in the IR region. As 
an appropriate object for comparison with our construction, we use the average 



A{Q) = ^ r dka^ik') (li; 
Q Jo 



Q 

that people manage to extract from jet physics data. Empirically, it has been claimed 
that this integral at Q — 1 -i- 2 GeV turns out to be a fit-invariant quantity. For it there 
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is an estimate: A{2 GeV) = 0.52 ±0.10 Our results for A{2 GeV) obtained by the 
substitution a^J^ and a'^ into Eq. (|ri|) for some values of the running coupling at the 
normalization point /i^ = are summarized in the Table. 



aan(M2) 


0.34 


0.36 


0.38 


Ai_ioop(2GeV) 


0.50 


0.52 


0.55 


A2_ioop(2GeV) 


0.48 


0.50 


0.52 



Note here that a nonperturbative contribution, like the second term on the l.h.s. of 
Eq. (H), reveals itself even at moderate Q values by "slowing down" the rate of as{Q) 
evolution. For instance, in the vicinity of the c and b quark thresholds at Q = 3 GeV 
it contributes about 4%, which could be essential for the resolution of the "discrepancy" 
between "low-Q" data and direct Zq measurement for as{Mz)- 

In this letter, we have argued that a possible way to resolve the ghost-pole problem for 
the QCD running coupling can be found by imposing the Kallen-Lehmann Q^-analyticity 
which reflects the causality principle of QFT. The analytic behavior in the IR region is 
restored by a nonperturbative contribution. The procedure of constructing the analytic 



running coupling is not unambiguous [jT2|. We have considered the simplest way which 
does not require any additional parameters and operates only with A or a value of the 
coupling at a certain normalization point. The requirement of analyticity yields significant 
modification in the IR and intermediate domains and leads to the universal value of 
Q;an(0). In this paper we have obtained the analytically-improved result aan(Q^) for the 
QCD running coupling that turns out to be quite stable with respect to higher-order 
corrections for the whole interval of and agrees with low energy experimental evidence 
for the IR-finite behavior. 

Our construction does not contain adjustable parameters. This is due, in particular, to 
the convergence of a nonsubtracted spectral integral, that is, with the asymptotic freedom 
property. Here, analyticity plays the role of a bridge between regions of small and large 
momenta. The idea that "analyticity is the key factor relating high energy to low energy" 
has recently been emphasized by Nishijima |]l3l in the context of connection between the 
asymptotic freedom and color confinement. However, in our approach this connection is 
not so direct. If, e.g., we admit (see Ref. fl^) the possibility of a UV fixed point for the 
QCD effective coupling at some small value of as < 0.05, then we arrive at a modification 
of the UV behavior with a power instead of (lncr)~^ decrease of the spectral function. So, 
in that case we can also use a nonsubtracted spectral representation. 

As far as it is difficult to present an explicit analytic expression for a^^iQ^), for the 
need of QCD practitioners, we propose an approximate formula. It can be obtained by 
the method of subtraction of unphysical singularities if one takes into account the explicit 
expression for the term that removes the ghost pole and an approximate expression in the 
form of the first term of the power expansion for the term that removes the unphysical 
cut. The corresponding expression reads 



1 

+ 7^ 



1 



2 1 -gvA2 



:i2) 



where ORciQ'^) defined by (|^) and, for n/ = 3, Ci = 0.035. Expression (0) ap- 
proximates the two- loop analytic coupling with an accuracy less than 0.5% in the region 
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2.5 A < Q < 3.5 A, practically coincides with the exact formula for larger values of mo- 
menta, and, therefore, can be used in the analysis of many experimental data. If there 
is necessary to consider a very small momentum, like Eq. (0), we suggest another ap- 
proximate formula which can be written in the form of Eq. (|^) with substitution, instead 



Q^/A^, the expression exp[ln Q^/A^ + bi In yln^ Q^/A^ + Air"^]. For Q > A, the accuracy 

of this approximation is less than 5% and yields only 3% error into the A{2 GeV) value. 
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